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^^ ' Abstract We present results of a population synthesis study for semidetached short orbital period binaries which contain low- 

c7^ I mass ( i 1 .5 Mq) donors and massive (> 4 Mq) compact accretors, which in our model represent black holes. Evolution of these 

^O . binaries is determined by nuclear evolution of the donors and/or orbital angular momentum loss due to magnetic braking by the 

^p ' stellar wind of the donors and gravitational wave radiation. According to our model the estimated total number of this type of 

fH ^ black-hole binaries in the Galaxy is t 10000. If the magnetic braking removing angular momentum in semidetached black-hole 

O ,■ binaries is described by the Verbunt & Zwaan formula, the model predicts that among them around 3000 systems with periods 

> 2 hours would be transient. In addition one finds several hundreds of luminous and stable systems with periods between 3 and 

8 hours. Several dozens of these bright systems should be observed above the RXTEASM sensitivity limit. The absence of such 

systems implies that upon Roche-lobe overflow by the low-mass donor angular momentum losses are reduced by a factor more 

— j I than 2 with respect to the Verbunt & Zwaan prescription. In such a case short period black-hole soft X-ray transients (SXT) 

• • ■ may have donors that overflow the Roche lobe in the core hydrogen-burning stage. We show that it is unlikely that the transient 

^ j^ ' behaviour of black-hole short-period X-ray binaries is explained by the evolved nature of the stellar companion. On the other 

k>( , hand a substantial fraction of black-hole binaries with periods > 3 hours could be faint with truncated, stable cold accretion 

> I ■ discs as proposed by Menou et al. Most of the semidetached black-hole binaries are expected to have periods shorter than ~ 2 

C^ ' hours. Properties of such, still to be observed, very small mass-ratio iq = M2/M1 < 0.02) binaries are different from those of 

their longer period cousins. 
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1 . Introduction ll998tlErgma & Sarnal200ll) . It is therefore of interest to check 

which models of binary evolution are capable of reproducing 

All known low-mass X-ray binaries containme black hole com- ,, , ... j, ., t iv/miin 1 .• j tu • u 

-' ° the transient nature or the LMBHB population and their ob- 

ponents (massive compact accretors, hereafter LMBHB) are j/jjj^ u .■^ . ^ .t 

^ ^j ^^ II ' I ' served (or deduced from observations) mass-transfer rates. In 

transient JMcClintock & Re millard 2005). The instability re- . ,.,. ... ., . . , , j . . ^ .u 

' ■' addition population synthesis models can be used to test the 

sponsible for this, often intermittent, behaviour is apparently ^ .• .• ^ .u j 1 ^u i. .• j .u 

^ I — rr .. j statistics of the model versus the observations and the existence 

the same that drives dwarf-nova outbursts (see lLasotall2001L f , . ., , n » ui 1 .• ^t a/tdod » 1 . ju 

^ ' ^' • of a hidden (quasi)s table population of LMBHB postulated by 

for a detailed review of the "disc instability model" - DIM) fTf r~j1 ITqqnl 

but in the case of LMBHB UTadia t ion of the outer disc play s 

a crucial role JvanParadiislll996t JPubus et al.lll999l EoOl]) . ^"^""g *e 17 confirmed black hole binaries, (i.e. 

For this instability to operate, rather low mass-transfer rates ^"*^ ^ measured mass-function larger th an ~ 3 Mn) ten 

are required ( < 10'° - lO" M^ yr'). This conjecture has are short-period (P„^ < Iday) systems j Orosz et alj |2004t 

led to the suggestion that all LMBHB are formed via a pe- NcClintock & Remillar j |2005|). Their evolution should be 

,. , .. ^u ■ u- u ^u A t » » similar to that of cataclysmic variables (CV) which prompts to 

culiar evolutionary path, in which the donor stars start mass -' \ / f f 

-,..,,,,,, 1 » .u A f ^u ■ ■ apply to LMBHB the model exploited for CV: evolution driven 

transfer to the black ho le only at the very end of their main- ^^ ^ ^ 

,-f .■ f \v t 1 llif.n/cl Fc p~c~j I mainly by angular momentum loss via magnetically coupled 

sequence lifetime (e.g. IKing et al.l 119961 lErgma & Fedoroval .? .7 & & .7 f 

stellar wind (MSW) and gravitational waves radiation (GWR). 

Send offprint requests to: J.-P. Lasota, lasota@iap.fr This is a non-trivial exercise because of the well-known prob- 
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lems of the "standard model" of CV evolution, e.g. the diffi- 
culty in reproducing the period gap and period mini mum ob- 
served in the distribution of CVs (see e.g. lKolbl2002i) . 

The criteria for the instability in the DIM depend mainly 
on the mass transfer rate, which is set by the period of the sys- 
tem and mass of the accreting body, i. e., along its evolutionary 
track a binary may move trough stable and unstable states. 

Therefore the evolution of LMBHBs provides an interest- 
ing challenge to our understanding of both the stellar evolution 
and the accretion-disc physics. This prompted us to carry out a 
population synthesis study for low-mass binaries with massive 
companions based on a systematic investigation of the evolu- 
tion of the latter systems under different assumptions on the an- 
gular momentum loss from the system. As the next step we at- 
tempted to infer the behaviour of a constructed population with 
respect to the stability criteria in DIM and to verify, whether it 
is possible that the Galaxy harbours a population of LMBHB 
with faint cold and stable discs. Population synthesis results are 
presented in Section|2and compared to observations in Section 
13 Discussion and Conclusions follow. 



2. Formation of low-mass black-hole binaries and 
their current Galactic population 

2.1. Formation of binaries containing a black hole and 
a low-mass star 

The formation of close binaries containing a black hole and 
a low-mass star is is still not well understood. The progeni- 
tor of the black hole must have been a very massive star, with 
heavy mass loss, which went through a supernova explosion. It 
is remarkable that apparently a low-mass companion can sur- 
vive such violent evolution. The "standard" formation chan- 
nel for such binaries is an extension of the original formation 
scenario pro posed for low-mass X- ray binaries with a neutron 
star accretor dvan den Heuvelll983l) : a massive star with a low- 
mass companion evolves and the low-mass star spirals into the 
giant in a common-envelope phase. The product of common- 
envelope stage is a compact binary consisting of the helium 
core of the massive star (a Wolf-Rayet star) and the low-mass 
companion. Neither the small amount of matter lost in super- 
nova explosion in which the neutron star is formed nor a na- 
tal kick imparted to neutron star disrupt the system, so the bi- 
nary survives and, due to angular momentum losses, becomes 
tighter and tighter until mass transfer ensues. In order to end 
with a black h ole, one simply needs a more massive star to 
start with (e.g. ' McClintock & Remil lard"l986!: Ide Kool et all 
Il987l) . Unfortunately, the above scenario contains many very 
uncertain aspects: the amount of mass loss of massive stars, 
the outcome of the common-envelope phase, the amount of 
mass loss of the helium core after the common-envelope phase 

and the mass loss and possible asymmetric k i ck induced at the 
^ I ' II I I II I 

supernova explosion (see e.g. jL anger 1989; Ibenetal. 1995; 
Wellstein & LangeJll999t lKalogera.1999: .Fryer & Kalogera 



2.2. Evolution of black-hole low-mass X-ray binaries: 
Evolutionary computations 

In order to be able later to compare the properties of the ob- 
served black-hole low-mass X-ray binaries with our model, we 
have made detailed calculations of a grid of evolutionary se- 
quences for different black hole and companion star masses 
and different initial periods. Depending on these parameters, 
the low-mass companions can overflow (if ever) their criti- 
cal lobes at different evolutionary stages. Consideration of a 
range of initial periods is important, in particular to assess 
suggestion than transient low-mass black-hole X-ray binaries 
must all harbour d onor-stars that almost exhausted hydrogen 
in their cores, e.g. (" King et al]ll996HErgma & Fedoroval 19981; 
Ergma & Sarna 200]^ 

Evolutionary confutations were carried out by means of 
the TWIN version of Eggleton (1971, priv. comm. 2003) evo- 
lutionary code. The most important updates of the code con- 
cerning equation of state, nuclear rea ction network and opac- 
ity are described by Pols et alJ (1199 5). Our version of the code 
implements angular momentum loss (AML) via a magnetically 
coupled stellar wind (MS W), following the model suggested by 
IVerbunt & Zwaan (^1981) and via gravitational wave radiation 
(GWR, Landau & Lifshitz 1971) . We should mention specifi- 
cally that, at difference with other codes which include AML 
via MSW, it is assumed that a critical minimum depth of the 
convective zone Z^om > 0.065 R© is required for sustaining a 
significant surface magnetic field. Within assumptions used in 
the code, this limits the masses of main-sequence stars which 
may have MSW to lower than 1 .6 Mq. On the other hand, since 
it is believed that the magnetic field is anchored in the radia- 
tive core, MSW is switched-ofF as soon as convective envelope 
penetrates to 0.2/?star ("disrupted" magnetic braking model). 
For initially unevolved MS-stars MSW terminates when their 
mass decreases to ~ 0.3 Mq. 

The Verbunt & Zwaan (19811) AML law can be written as 



dJ_ 

dt 



-Q.5lQ'^^r^k^MdRW dyncm, 



(1) 



200lUNelemans & van den Heuyei200lHjustham et al.l2006l) . 



where k^ ~ 0.1 is the gyration radius of the donor star and M^, 
Rd and o) are respectively its mass, radius and rotation angular 
frequency, / is a parameter calibrating the strength of MSW 
AML. In the evolutionary code, the gyration radius is estimated 
for every model. 

Solar chemical composition of the models was assumed. 
Trial computations for donors with metallicity Z=0.008 did 
not reveal any serious dependence of results on Z. Different 
metallicities might also influence the mass of the black hole, 
as the mass loss rates of massive stars are expected to scale 
with Z^-^. However, what really matters for our purpose is the 
pre-SN mass. Since both mass loss rates for MS- and WR- 
stars are highly uncertain we neglect the Z-dependence of stel- 
lar winds. Comparison of our pre-SN masses wit h the results 
of trial computations with the iHurlev et alJ ( 120001) code which 
employs a different combination of stellar winds shows that 
pre-SN masses are the same within a few solar masses for 
metallicities of 0.02 and 0.008. We neglected the possi ble influ- 
ence of an irradiation induced stellar wind ( Tavani & London! 



L.R. Yungelson et al.: Low-mass binaries and SXT 



19931) on the donor evolution. This is justified since irradiation 
feedback is important o nly for systems with M ~ Mq donors 
llBuning & Rittenl2005l) which comprise a negligible minority 
of population under consideration. 

In codes describing CV evolution the AML is usually nor- 
malized in a way to obtain the correct location of the period 
gap (see e.g. Kolb & Willems 2005) . Since no gap has been 
detected for the LMXB orbital distribution we left our code 
uncalibrated. Let us add that at variance with cataclysmic vari- 
ables, in systems with massive primaries the donors are only 
slightly out of thermal equilibrium and the width of the calcu- 
lated "period gap" is small, compared to CVs. For the models 
that overflow Roche lobes when they have already burnt a sig- 
nificant fraction of hydrogen in their cores the gap disappears 
(see Fig.0]below). 

The grids of evolutionary models were computed for ini- 
tial secondary masses (M20) of 0.4 - I.6M0 and initial black 
hole masses (Mio) of 4, 7, and I2M0. In the first set of com- 
putations a "standard" model of evolution of low-mass binaries 
with AML via MSW (/ = 1 and disrupted magnetic braking) 
and GWR was assumed (henceforth, model A). Since popu- 
lation synthesis under these "standard" assumptions did not 
provide satisfactory results (see below), in the second set of 
computations we assumed ad hoc that MSW does not operate 
in semidetached systems with black-hole components (hence- 
forth, model B). In model B only GWR operates after RLOF. 
We also used for comparison a model similar to model A in 
which the MSW angular-momentum loss rate operated per- 
manently, but reduced by a factor ~ 2 if - 1.37 - model 
C). Evolution models A, B, and C give rise to 13737, 12089, 
and 11659 short-period semidetached systems, respectively. 
Corresponding LMBHB -population models are compared in 
detail in Sec.|3 

The models for low-mass donors typically ceased to con- 
verge at M™" ^ 0.03 M0 in evolutionary tracks for models A 
and C and at M™« ^ 0.06 Mq for model B. This happens be- 
cause of the absence in the code of adequate tables of EOS 
and opacities for that low-mass stars the slight difference in the 
limiting masses stems from some differences in the versions of 
the code used for models A and B, C; in our results this only 
slightly influenced M™" and position of orbital period mini- 
mum). For population modeling, all tracks were truncated at 
M2""=O.O2M0 or 0.06 Mq, respectively. However, already at 
M2=O.O6M0 the mass ratio of components q < 0.02, the cir- 
cularization radius of accretion stream for such low q becomes 
greater than the outer radius of the accretion disc and it is un- 
clear, whether mass transfer happens at all (see 13.31 for more 
details). Thus, although evolutionary computations were con- 
tinued beyond q ^ 0.02 with the same formalism as for larger 
q, the results may be not applicable to any physical systems. 

2.3. The Galactic population of low-mass binaries with 
black-hole companions 

We modeled the Galactic population of black-hole 
binaries using the population synthesis code SeBa 
llPortegies Zwart & Verbunn 1 19961) with updates described 
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Figure 1. Initial mass - final mass relations for massive stars 
in binaries. The lines show the mass of stars at the end of the 
main sequence (TAMS), the initial mass of helium core Mcore.i, 
the pre-supernova mass Mpre-SN, and the final black-hole mass 
Mbh- 



in IPortegies Zwart & YungelsonI ( Il998h : iNelemans et alJ 
('20011): iNelemans et al.l (f2004). The important points to note 
are that we assum e a 50% binar ity rate (2/3 of stars in bina- 
ries), an IMF after iKro upa et alJ lll993). an initial distribution 
of separations in binary systems (a) flat in logo between 
contact and 10*' Rq, a flat mass ratio distribution, and initial 
distribution of eccentricities of orbits H(e) = 2e. 

A summary of the assumptions about the evolution of mas- 
sive stars and Wolf-Rayet stars that we use in our code is shown 
in Fig. n where, for components of close binaries with initial 
masses of Mq = 25-100 Mq, we show their masses at the end 
of the main sequence (TAMS), their initial helium core mass, 
their pre-supernova mass and the final black-hole mass'. We 
assume that the stellar wind mass-loss rate increases in time 
and that in their total lifetime stars may lose an amount of mat- 
ter that increases with initial mass: O.OIM^. If the star loses its 
hydrogen envelope, we switch to the Wolf-Rayet stellar mass- 
loss prescription. Mass loss on the main-sequence causes stars 
with masses of about 85Mq to lose all their hydrogen before 
they reach the TAMS. For more massive stars the mass loss 
is so uncertain that we assume for simplicity that they lose 
42 Mq on the main sequence, leading to the peculiar upturn 
of the solid line in Fig. [2 However, these extremely massive 
stars are so rare that they do not contribute to the population of 
low-mass black-hole binaries. Fo r the common-envelope phase 
we us e the standard prescription (Webbink 1984:' de Kool et alJ 
Il987l) . with the efficiency and structure parameters a and A 
combined: aA = 2. For the mass loss by Wolf-Rayet stars 
formed af ter the common-envelope phase, w e use the law de- 
rived by .Nelemans & van den Heuvell ( 1200 ih . which is based 



' The overwhelming majority of close binaries donor stars lose their 
hydrogen envelope through mass transfer after the main sequence, but 
before the onset of helium core burning. For such binaries the final 
masses are rather insensitive to the initial period of the binary and the 
initial q. 
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on the compilation of estimated mass-loss rates of observed 
Wolf-Ray et stars of lNugis & Lamerj ilOOd) . 

The formation o f black holes in the code follows the fall- 
back prescription of iFrver & Kalogeral 11200 11) whose basic as- 
sumption is that a constant fraction of the supernova explosion 
en ergy is used to u nbind the envelope of the star At difference 
to iFrver & Kalogera C2001.) . who consider the explosion en- 
ergy as a function of pre-supernova mass, we keep it fixed at 
10^" ergs, which is within the expected range but favours for- 
mation of rather massive black holes (up to ~ 15 Mq). For sys- 
tems that are relevant for this study, i.e. naked helium cores, the 
prescription we use actually almost invariably results in com- 
plete collapse of the pre-supernova object: the CO core is very 
tightly bound while around this core only a small layer of he- 
lium is left because of the strong stellar wind during the pre- 
ceding Wolf-Rayet phase. We assume tha t black holes receive 
a small asymmetric kick at formation (e.g. | jonker & NelemansI 
12004'; Willems et al. '2005"; Gualandris et al. 2005), which is 
taken from a Paczvhski (1990) and Hartman (1997) velocity 
distribution with a dispersion of 300 km s~\ but scaled down 
with the ratio of the black hole mass to a neutron star mass. 
The population of short-period low-mass black-hole binaries 
is not sensitive to the assumed kick distribution, since scaled 
down with Mbh/Mns -ratio kick velocities are too small to dis- 
rupt close binaries in SN explosions. This was confirmed by 
a test run assuming Maxwellian kick d istribution for pulsars 
with cr - 265km/s (Hobbs et al. 2005). The total number of 
systems that filled Roche lobe in Hubble time and evolved to 
shorter periods (in models A and B, see below) changed by 
several dozens only. 

The model of the Galactic population of detached low-mass 
binaries with black-hole companions is then produced by fold- 
ing the results of approximate evolutionary computations for a 
set of 250000 initial binaries with Mio > 25 M© with a time 
and position dependent model for the st ar formation history 
of the Galaxy that is based on results of iBoissier & PrantzosI 
( 1199 9') and a model for th e distribution of stars in the Galaxy 
(see Nelema ns et al.ll2004l) . Some details of this model are de- 
scribed in AppendixlAl 

With our assumptions we find that, within a Hubble time 
(13.5 Gyr), ~ 49 000 binaries that have orbital periods be- 
low 2.0 day and contain black holes accompanied by main- 
sequence stars less massive than 1.6 M© were formed in the 
Galaxy.^ We restrict ourselves to M20 < 1 .6 Mq since for more 
massive stars magnetic braking is not supposed to operate. Out 
of these ~ 49 000 binaries, ~ 17,000 were brought into con- 
tact in a Hubble time, but only ~ 14000 evolved to shorter 
periods (Table O. If only AML via GWR acts in systems with 
more massive donors, they experience case AB mass-transfer-', 
evolve into subgiants with increase of orbital periods and end 
their lives as helium white dwarfs. 




5x10 



t, yr 

Figure 2. Solid histogram - fractions of low-mass bh-i-ms bi- 
naries formed at different epochs in the lifetime of the Galaxy. 
Dotted curve shows the shape of star formation history adopted 
in the population synthesis code. At maximum, the star for- 
mation rate (SFR) is ^ 15.6 Mgyr '. Current SFR is ^ 
1.9 Moyr-i. 



Because black holes form in the first ~3 Myr after their 
progenitor formation, the formation history of bh-i-ms binaries 
(Fig.O strictly follows the star formation history. 

The upper panel of Fig.|3lshows the relation between initial 
masses of MS -components and orbital periods after circulariza- 
tion of the orbits after black hole birth events'*. The lower panel 
represents relation between masses of components in these sys- 
tems. 

2.4. The current Galactic population of black-hole 
low-mass X-ray binaries 

We have simulated the current Galactic population of semide- 
tached low-mass black-hole binaries by convolving the above- 
described "underlying" population of binaries born at differ- 
ent epochs in the history of the Galaxy with the grids of pre- 
computed evolutionary tracks for low-mass components in bi- 
naries with different combinations of components and post- 
circularization (initial) orbital periods. 

Since the features of the model population may be under- 
stood as a consequence of different patterns of evolution of low- 
mass donors (Fig.|4} we will first present examples of evolu- 
tionary tracks. The parameters of the tracks shown in Fig.|3]are 
hsted in Table [2 For every combination of M20 and M^q there 
are systems that, upon Roche lobe overflow (RLOF), evolve 
to shorter periods having mass-transfer rates M x (10^^ - 
10"'") Moyr"'; for M20 - i- Mq this happens for systems 



^ This number is ~0.01% of all binaries with black-hole compo- 
nents formed in Hubble time. 

' According to common classification, in "case AB" mass-transfer 
starts when the donor is still in the core-hydrogen burning stage and 
continues in the hydrogen- shell burning stage. 



'^ Because of specifics of time-realizations in the population synthe- 
sis code every dot in the plot actually represents four or five systems 
with similar combination of masses and orbital period, but bom at dif- 
ferent time. Since we present one random realization of the model for 
a population of black-hole+main-sequence star systems, all numbers 
given are subject to Poisson noise. 
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Figures. Upper panel: relation between orbital periods and 
masses of main-sequence companions of black holes after cir- 
cularization of the orbits. Systems with main-sequence com- 
ponents that may fill their Roche lobes in their main-sequence 
lifetime or in the Galactic disc lifetime (13.5 Gyr) are located 
to the left of the solid curves. Left and right solid curves cor- 
respond to Mbh - 4 Mq and 12 M©. Lower panel: relation be- 
tween initial masses of components in bh+ms systems. 



with an abundance of hydrogen in their center Xc ^ 0.05 at 
RLOF. These systems may have post-circularization periods 
up to ~ 40 hr (Fig.O and they overflow Roche lobes at pe- 
riods i 12 hr (see Tableland Fig. |4} since the period of 
semidetached systems depends only on the mass and radius of 
the donor and the latter changes only slightly in the hydrogen 
burning stage or practically does not change in Hubble time for 
M20 i 0.95 Mq. Systems with hydrogen-exhausted cores (two 
rightmost tracks in Fig. 01 No. 7 & 8 in Table[Q evolve to longer 

S ieriods. This is the well known effect of "bifur cation period" 
Tutukoveta l."l985': 'P vlvser & Savoniie"l989'). The bifurca- 
tion period Pb defines the right border of the position of new- 
born systems that evolve to shorter orbital periods upon RLOF, 
as shown in Fig.|2 

In systems that are initially driven by both MS W and GWR, 
M drops rapidly from M - 10"' Mq yr"' to ~ 10"'° Mq yr"' 
after termination of MSW. Like in cataclysmic variables, the 
evolution continues with decreasing M and there exists a min- 
imum orbital period ~ (70 - 80) min. LMBHB close to the 
period minimum have yet to be observed but our population 
synthesis model predicts that they might form the majority of 
LMBHB. These systems have very low mass ratios q < 0.02 



Table 1. Parameters of tracks shown in Fig. ^ The columns 
list initial orbital period of the system, period at RLOF, central 
hydrogen abundance at RLOF, mass of helium core at RLOF. 
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Figure 4. Mass transfer rates vs. orbital period for systems 
with initial masses of components 1 M© and 7 M© that have 
different periods Po after circularization of orbits. Solid lines 
show tracks of donors evolving under "standard" assumption 
that AML via MSW operates until secondaries become con- 
vective (/=1). The numbers along the tracks correspond to 
the numbers in Tabled The leftmost line is for Pq = 9.6 hr; 
in this system secondary fills its Roche almost unevolved at 
P ^ 7.44 hr Note the presence of a tiny period gap for this 
system. The systems to the right have Pq from 38.4 to 45.6 hr 
(respectively, at RLOF they have, either X^ < 7 ■ 10"^ or he- 
lium core with Mue < O.OIMq, see Table Q. The dotted line 
shows mass transfer rate critical for disc instability (Eq. |2j- 
The dashed line marks the approximate upper limit to the av- 
erage mass transfer rates for short-period black-hole binaries, 
as given by Eq. (|8}. The diamonds show estimates of mass- 
transfer rates based on assumed recurrence period, see Table 

m 



and it is far from obvious that mass transfer occurs in them in a 
similar fashion to standard compact semidetached binaries (see 
SectionlO. 

For the most massive stars that still have MSW and fill their 
Roche lobe close to Pi„ nuclear evolution continues to play a 
certain role after RLOF. As a result, they first evolve to shorter 
periods, but, when hydrogen is almost completely exhausted 
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in their cores, turn to longer periods (tracks with "loops" in 
Fig. 0}. In the extreme case of M20 = 1.6Mq donor, this "u- 
turn" may happen at a period as short as ~ 4 hr, but such sys- 
tems with periods from a very narrow interval of initial Porb are 
extremely rare in the parental population. 



2.4.1 . Evolved companions ? 



The systems that have at RLOF 0.05 i X^. i evolve to short 
periods and may have from the very beginning of the Roche- 
lobe overflow M < 10"' Mq yr"'. In the period range 4 - 
12 hr they have rather low mass-transfer rates: M ~ 10""^ - 
10"" Mq yr"' ( iTutukov et al. 1985). They are represented by 
the tracks No. 3-6 (with initial periods 40.8 - 43.2 hr and pe- 
riods at contact 12.5 - 14.8 hr) in the centre of Fig.|3 In the 
model population, they do not evolve to periods below about 4 
hr, since they have to be already significantly evolved at the in- 
stant of RLOF which leaves them only a few Gyr of evolution- 
ary lifetime in semidetached state. Since the interval of periods 
of these systems and their mass-transfer rate M overlap with 
the those of SXT it is tempting to assume that they provide 
the explanation of the properties of these systems. However, 
this cannot be the case because in the majority of systems the 
donors fill their Roche lobes almost unevolved: for the systems 
shown in Fig.|3 the ratio of initial period intervals of systems 
that fill Roche lobes "evolved" and "unevolved" is ~ 8%. The 
actual ratio of the numbers of stars in two groups is much lower, 
since only systems with M20 i 0.95 M© may exhaust hydro- 
gen in their centres in Hubble time. This is confirmed by the 
upper panel of Fig.|5l(see below), which shows that even the 
upper estimates of the average mass-transfer rates in observed 
SXT hardly may be explained by the standard AML via MSW 
and GR model, because of low number of systems with desired 
initial parameters in the parental population and their short life- 
time in the observed period range. Simply removing all systems 
that will fall in the "unevolved" class by fine-tuning the initial 
period distribution of bh+ms is difficult: the most obvious way 
would be to change the common-envelope efficiency. However 
we tried this by repeating our calculations with aA - 0.5 and 5 
(the standard value in our model is 2), but the "unevolved" class 
still dominates. The number of black hole + low-mass MS-star 
systems formed in a Hubble time was, compared to "standard" 
model, 50% lower and 15% higher respectively for these two 
calculations.^ 



This shows that the transient nature of the short pe- 
riod LMBHB most likely cannot be explained by the pres- 
ence of stron gly evolved main-sequence com panions as pro- 
posed by e.g. iKing et alJ ( 1 1996.) : .Ergma & Fedoroval (Il998.) : 
lErgma&Sarnal ( l200lh . 
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Figures. Mass transfer rates vs. orbital period for the model 
population of semidetached low-mass binaries with assumed 
black hole accretors. Upper panel shows systems that have ir- 
radiated discs that are stable according to criterion given by 
Eq. (|2jl in the model for a "standard" assumption on AML 
via MSW always operating in stars with convective envelopes 
(/ = 1, Model A). Middle panel shows population of systems 
with unstable irradiated disks in the same model. Lower panel 
presents population for the model in which MSW does not op- 
erate in semidetached systems with black hole accretors (model 
B). In this model all binaries have unstable irradiated disks. The 
break-down over number of stable and unstable systems in par- 
ticular models is listed in Table |2l The dashed line marks the 
approximate upper limit to the average mass-transfer rates for 
short-period black-hole binaries, as given by Eq. ^, while di- 
amonds show estimates of rates based on assumed recurrence 
period of 30 years, except for V616 Mon and IL Lup whose re- 
currence times are known to be ~ 50 and 10 years, respectively 
(see Table|3]and text for details). 



3. Accretion in tlie modeled LIVIBHB population 
and its observational consequences 

3.1. Bright, persistent systems 

3.1.1. Disc stability 

Figure |5] (upper and middle panels) shows the modeled popu- 
lation of semidetached low-mass binaries with assumed black 
hole accretors in "standard" model (A) with disrupted mag- 
netic braking. In Figure [Sj only systems for which DIM is 
valid and stability analysis may be applied are plotted (~ 3000 
systems). The remaining systems (Table |2j have low-mass ra- 



Respectively 24,000 and 58,000 systems. 
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tios of components (q < 0.02) and low-mass (M2 < 0.03 Mq) 
donors. If mass transfer occurs and can be described by the 
"standard" model of evolution driven by AML by GWR, these 
systems should be concentrated around a "period minimum" 
at 70 - 80 min and beyond it and at mass-transfer rates 
10"'° - 10"" Moyr"'. Model A contains ~ 13000 systems 
(Table ^ of which some may be luminous and non-transient 
(upper panel of Fig.|5jl. Such steady systems should satisfy the 
stability criterion according to which irradiated discs are sta- 
ble if their accretion rate exceeds certain critical value M^.^;, 
JvanParadiisll 19961) . According to lPubusetTal] ( 1 19991) the lat- 
ter is 



M., 



2.4x10""Mb2['^ 



1 1010 cm/ 



C 



5 X 10-4 



-0.5 



Mo yr" 



Table 2. Break-down of black-hole candidate populations into 
sub-groups according to accretion disc stability. Model A - 
"standard" evolution model with disrupted magnetic braking 
(/ = 1), Model B - model without MSW in semidetached sys- 
tems. Model C corresponds to AML by MSW reduced by a 
factor 2 (/ = 1.37). The total number of systems in contact is 
indicated for each model. The population is then separated into 
systems with q < 0.02 for which description in all evolutionary 
calculations may be inadequate (see Section 031 and systems 

, with q > 0.02. For the latter, the systems are further broken 
down into bright persistent systems (DIM stable, Section ITTl 
and transient systems. Some of the systems marked a "tran- 
sient" could be cold, stable, truncated-disc systems (MNL sta- 

(2)ble, Section lTT^ . 



where Mbh is the mass of accretor in Mq, Rij is the disc outer 
radius, and C is a measur e of heating of the disc by X-rays 
llShakura & Sunvaevlll973l) : for a given disc radius R the irra- 
diation temperature ![„ is 



C 



Atto-R^ ' 



(3) 



where cr is the Stefan-Boltzraann c onstant. We use C = 5 x 
lO^i, following IPubus et alJ 1119991) who have found that this 
value is consistent with properties of persistent low-mass X- 
ray sources. One should, however, keep in mind t hat C does 
not ha ve to be a constant and might e.g. vary in time (lEsin et alJ 

knod . 

The disc outer radius Rj, which depends only on the mass 
ratio once scaled to the or bital separation, is estimated using 
Tab. 2 of IPaczviiskjl ( 1 19771) who computed the largest possible 
stable, periodic non-intersect ing orbits. This estimate is v ery 
close to the one obtained by IPapaloizou & Pring ld (119771) by 
studying the linear response to tidal forcing. For mass ratios 
q X 0.02, the disc outer radius becomes comparable to the 
circularization radius Rdic and the Rd estimate is not valid. It 
is then uncertain how accretion proceeds. In what follows we 
concentrate on systems with mass ratios q > 0.02 (3272 sys- 
tems in model A), corresponding to periods longer than ^ 2 
hours, and come back to the systems with low mass ratios in 
SectionlO 



3.1 .2. Excess of luminous persistent LIVIBHB 

The number of systems for each model population is given in 
Table |2 which also gives the number of luminous, persistent 
systems expected in each case according to the criterion 0. 

The consistency of these numbers can be compared to ob- 
servations, which show no persistent LMBHB with short Porb. 
Possible exceptions, i.e. systems showing variability that could 
be different from what is called a transient behaviour, such as 
GX 339-4 or GRS 1915-H105 have long orbital periods and 
are not addressed here. Model A (standard MSW) predicts the 
highest number of persistent sources. These 367 systems have 
an average luminosity L - O.lMc^ ^ 5 • 10^^ erg s'. Given 
their location in the Galaxy provided by the population synthe- 
sis code, we estimate the intervening absorption assuming that 
the gas is, roughly, distributed as « = 1.8exp(-t//3.5 kpc - 





Model A 


Model B 


Model C 


Total 
q > 0.02 


13737 
3272 


12089 
5074 


11659 
2991 


Hot stable 
Transient 
(Cold stable) 


367 

2905 

(62) 




5074 
(2495) 


107 

2884 

(79) 



z/0.125 pc) cm"^ with d the distance to the Galactic centre, z 
the height above Galactic plane and where the gas is distributed 
like the thin disk population (Bahcall 1986) with a density and 
scale-height taken from Zomb eckl ( 11990 ). The average column 
density along the line of sight from the Sun (at 8.5 kpc from 
the centre) to the systems is logNn » 22.1. 

The accretion luminosity may be assumed to be emitted as 
that of a 1 keV blackbody, corresponding to the 'high/soft' X- 
ray spectral state of LMBHB. Taking into account X-ray ab- 
sorption, the number of persistent systems expected above a 
threshold of 7 -lO"" erg s"' cm"^ in the 2-10 keV band is 
85. This threshold corresponds to the completeness limit for 
the X-ray sources seen by the RXTE All Sky Monitor (about 3 
mCrab or 0.2 counts/s) deduced from the log N - log S distri- 
bution of extragalactic sources bv lGrimm et alJ ( 120021) . Setting 
the threshold ten times higher still yields 13 bright, persistent 
LMBHB with ASM fluxes comparable to e.g. SMC X-1. 

Assuming a power-law spectrum with photon index -2 (flat 
in vFy), corresponding to the 'low/hard' X-ray spectral state 
of LMBHB, one obtains almost the same number of detectable 
persistent systems: a 1 keV blackbody peaks at 2.7 keV so the 
power-law emission is not much different. Persistent LMBHB 
might escape detection if they emit most of their energy below 
1 keV. Indeed, KV UMa in the low/hard state had blackbody- 
like emission peaking at ^ 25 eV, with a spectral luminos- 
ity almost 10 times greater than in X-rays wh ere it displays a 
power law spectrum (McClintock et al.ll200ll) . Assuming that 
only 10% of the accretion luminosity is emitted in hard X- 
rays would still imply 16 persistent systems above 3 mCrab 
in the RXTE ASM (1 above 30 mCrab). Alternatively, the discs 
may be not irradiated, in which case no systems would be sta- 
ble; this, however, would be surprising as both persistent NS 
LMXBSs and outbursting black hole transients show clear sig- 
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natures of irradiation. Note that with parameter C - 5 x 10 "* 
(Eq.|3l, only a very moderate 0.5% of the accretion luminosity 
(efficiency of 10%) is reprocessed in the disc. 

Despite obvious uncertainties, the large number of systems 
at high mass-transfer rates in Model A makes it difficult to 
avoid the presence of a significant number of observed per- 
sistent systems. However, mass-transfer rates could be smaller 
than suggested by models of binary evolution driven by AML 
through MSW and GWR. In such a case Model B (GWR only) 
predicts no persistent systems at all while the intermediate case 
represented by Model C (reduced MSW) yields 29 systems 
above 3 mCrab (or 5 above 30 mCrab) in the RXTE ASM. In 
addition mass-transfer rates predicted by model A are too large 
compared to observations of transient systems. These points are 
further discussed in the next section. 

3.1.3. Reduced Magnetic Braking 

While energy and angular momentum loss via gravitational 
wave radiation is beyond doubt, the issue of angular momen- 
tum loss by binaries via magnetically coupled stellar wind and 
its quantitative description are much more arguable. The earli- 
est formulation of Eq. Q by Verbunt & Zwaan ( 1981) and its 
variation s are based on ex trapolation of stellar rotation brak- 
ing law JSkumanichll 19721) over several orders of magnitude 
in stellar rotation rates, imply efficient spin-orbit coupling for 
low-mass star and has been modified and/or challenged by 
many authors Mameurv et al. 1988~ Collier CameronI l2002t 
Andronov et al. 2003 ; Ivanova&Taam 2003, and others); 
Baraffe & Kolbl (|2000J) argued that in order to explain the cata- 
clysmic binary period distribution (i.e. their evolution) a mag- 
netic braking law according to which the mass-transfer rate 
would, contrary to the "standard" case, decrease with increas- 
ing orbital period might be necessary. It is therefore fair to con- 
clude that the so-called "magnetic braking mechanism" is so 
uncertain that it can be almost modified at will. 

For this reason we computed a grid of evolutionary se- 
quences for binaries with low-mass donors and massive accre- 
tors under the assumption that AML via MSW does not operate 
in semidetached bh +ms systems. T his choice was motivated 
by stu dies such as lAndronov et alJ (i2003) : Ivanova & Taam 
1I2OO3I) according to which the MSW operating in close bina- 
ries is reduced by at least one order of magnitude with re- 
spect to Verbunt & Zwaan ( 198 1 ) law; but mass-transfer values 
obtained in this way should be rather treated as lower limits. 
Based on this set of evolutionary tracks we derived the second 
model for the population of semidetached binaries with black 
holes (model B). The model is presented in the lower panel of 
FigE] In this model the current Galactic population of semide- 
tached bh-Hms binaries contains in total, x 12000 objects, about 
5000 of which have q > 0.02 (TableEJ. 

The absence of MSW AML results in lower M in the 
systems with similar initial masses of components and post- 
circularization periods. As a result, in all systems of model 
B accretion discs are unstable according to the criterion (|2}. 
Instability limit is given by the lower bound of region filled by 
stable systems in "standard model". In a model without mag- 



netic braking in semidetached systems typical mass-transfer 
rates M are by factor up to 10 lower 

If MSW were to be absent from the very beginning, before 
the RLOF, the results would not change much. The number 
of systems would be reduced by a factor ~ 2. However, the 
longest period systems would be absent and this would make 
consistency with observations worse. 

The absence of stable bright low-mass sources argues in 
favour of an A ML reduced as compared to predictions of 
IVerbunt & Zwaaa (■1981.) magne tic braking model. It is worth 
noting that iMenou et al.l 1119991) also found that black-hole 
SXTs evolution should be driven by AML without MSW. In the 
view of both theoretical (the magnetic braking law, the value of 
C) and observational (the number of luminous, stable but unde- 
tected systems) uncertainties it would not make much sense to 
try to fix the value of parameter / in Eq. that would bring the 
calculated population in exact agreement with observations. An 
experiment with MSW angular momentum loss rate reduced 
by factor 2 (Model C, see Tabled) shows that finding such an 
agreement will require a more substantial reduction of the mag- 
netic braking AML. 

3.2. Transient systems and cold, stable systems 

The prediction of our model can also be compared to mass- 
transfer rates deduced from observations of SXT sources. Of 
course such deductions are notoriously uncertain but, as we ar- 
gue below, they can provide interesting information especially 
when combined with the knowledge that all observed short- 
period LMBHB are transient X-ray systems. 



3.2.1. Reservoir discs 

The usual method for SXT consists of obtaining the mass- 
transfer rate by dividing the m ass accreted during outburs t 
by the recurrence time (see e.g. IWhite & van Paradiislll99fl) . 
Regretfully, among short-period SXT the recurrence time is 
known only for A0620-00 (about 60 years) and 4U 1543-47 
(about 10 years). For the other systems one can only obtain up- 
per limits by assuming that since only one outburst has been 
seen the recurrence time is longer than 30 years. These esti- 
mates of M , derive d from the data on outburst parameters from 
IChen et alJ ill 997") and fo r KV UMa and V406 Vul based of 
observ ations presented in IChatv et al.l ll2003.) and iHvnes et alJ 
(120021) respectively, are given in the penultimate column of 
Tableland are shown in Figs.|3|5l and|6l 

Clearly even pure GWR cannot explain such low mass- 
transfer rates if one assumes that they represent secular val- 
ues. No LMBHB with a period i 4 hr can have a secular 
mass-transfer rate as low as M ~ 10"'^ Moyr"'; it would 
take post-minimum period systems several Hubble times to 
reach these periods and corresponding mass-transfer rates. The 
only known type of systems where such M may be expected, 
are binaries with initially degenerate dw arf donors, but eve n 
then orbital periods should be < 2 hr ( iDelove et al J 120051) . 
If the extremely low mass-transfer rate deduced for short pe- 
riod SXTs were due to downward fluctuations from an intrin- 
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Table 3. Orbital periods, mass-transfer rates based on recur- 
rence times and upper limits on the average mass-transfer rate 
for black-hole candidate short-period binaries. See text for de- 
tails. 



Name 



day Mo yr"' Mq yr 



.-1 

XTEJl 11 8+480 (KVUMa) 0.171 1.9 x 10"'^ 1.5 x 10^'" 

GROJ0422+32(V518Per) 0.212 1.3x10-" 2.3x10-'" 

GRS 1009-45 (MM Vel) 0.285 4.4 x 10"" 3.8 x lO"'" 

XTEJ1650-500 0.320 2.0x10"" 4.8x10-'" 

A0620-00 (V616 Mon) 0.323 3.3 x 10^" 4.8 x 10"'" 

GS 2000+25 (QZ Vul) 0.345 2.0 x 10^'" 5.4 x 10^'" 

XTE J1859+226 (V406 Vul) 0.383 4.1 x 10"'" 6.4 x 10"'" 

GRS 1 124-68 (GU Mus) 0.433 3.4 x 10"'" 8.0 x 10"'" 

H 1705-25 (V2107 Oph) 0.521 5.5 x 10"" 1.1 x 10^' 

4U 1543-47 (IL Lup) 1.125 4.2 x 10"'" 4.3 x 10^' 



sically high secular mass-transfer rate, one would still have to 
explain why the bright counterparts are not observed (Section 
13. U . Furthermore, some systems would then necessarily have 
higher rates than their secular value (by how much depends on 
the duty cycle of the purposed low states), hence there should 
be even more persistent systems. 

For longer periods, Fig.|4]shows, M of systems with Poib ~ 
8-10 hr may be consistent with evolution of systems where 
donors overflow Roche lobes extremely close to TAMS, but 
then one needs to assume some very special initial distribution 
of binaries over orbital separations which will finally, after sev- 
eral evolutionary stages with mass transfer and mass and mo- 
mentum loss from the system and supernova explosion, result 
in concentration of zero-age black-hole and main-sequence- 
star systems just at the desired very narrow range of separa- 
tions. 

The longest period SXT of the "short-period" group, 
4U1543-47 is, in principle, consistent with the evolution of 
systems with donors overflowing Roche-lobes after forma- 
tion of He-cores. But again, these binaries evolve through the 
Hertzsprung gap extremely fast, and the probability of observ- 
ing a descendant of such systems is very low, see discussion 
of intermed iate-mass "Hertzsprung gap" black-hole binaries in 
iKolhl ill 9981) . 

3.2.2. Truncated discs 

It is, however, plausible that the mass-transfer rates in short- 
period SXTs are much higher than those derived above. The 
method used to obtain the low values mentioned above treats 
the disc as a reservoir which during quiescence is filled up 
to a critical value at which the outburst is triggered and the 
disc emptied. The implicit assumption in this picture is that the 
reservoir is not leaky, i.e. that the accretion rate at the disc's 
inner edge is much smaller than the mass-transfer rate, as in- 
deed predicted by the DIM when the disc extends to the in- 
nermost stable circular orbit. However, there are many rea- 
sons to think that quiescent (and not only quiescent, see e.g. 



IPone & Gierliiiskil20Q6l) SXT discs are truncated and therefore 
"leaky" because of the inner trunca tion required to explain ob- 
served quiescent luminosities (e.g. lLasotalll99(Tl: iLasota et alJ 
1996). Also, Dubusetal. (2001) showed that the disc insta- 
bility model can reproduce the observed X-ray light-curves 
only if discs in SXTs are truncated (and irradiated). The trun- 
cated disc "paradigm" is now commonly used in describing 
the observed timing and spectra l properties of X-ray LMBHB 
( iMcClintock & Remillardl2005l and references therein). 

For systems with non-stationary quiescent accretion discs 
the mass-transfer rate can be estimated as ( lLasotal200ll) 

eA^D.max 



M, 



where 



Mn 



•Mi, 



= 2.7 x 10^'c 



-0.83 



Ml 



-0-38 ^^ 3.14 

1 10"' cm/ 



(4) 



(5) 



is the maximum quiescent-disc mass and e = AMo/MD.max the 
fraction of the disc's mass lost during outburst; a is the kine- 
matic viscosity parameter, Rd is the disc's outer radius and Min 
is the accretion rate at the disc's inner edge. The usual, non- 
leaky disc estimates, neglect Min. For a stable equilibrium disc 

frecc = oo and Mtr = Min. 

Ho wever, in SXTs Mn is in fact the dominant term as ar- 
gued by Menou et al. ( 1999). They applied models in which the 
inner part of the accretion flow was represente d by an ADAF 
JNaravan & Yilll994 lAbramowicz et aljri995l) to several qui- 
escent SXT systems. Fitting the model to the observed spectra 
provided estimates of Min = Madaf which are higher than the 
"refilling" term eMo.max/frecc, for which except for A0620-00 
(and the long period SXT system V404 Cyg) only upper limits 
co uld be obtained . 

iLasotalOOOOl) suggested a general estimate of Min in quies- 
cent SXTs. It is based on the fact that in quiescence according 
to the DIM the disc accretion rate has to satisfy the inequality: 



M{r) < M,rit ~ 6.4x10" 



Mo) 1 10"' cm/ 



M0yr-i(6) 



10/ \10"'cmy 

JDubusll 19981 iLasotalEoOll note that this criterion which re- 
quires the whole quiescent disc to be in a cold thermal equi- 
librium state is different from the instability condition of Eq. 
13. Therefore, the mass accretion rate at the truncation radius 
(Mn) must be smaller than M^^.i^(rin)- Inside the truncation ra- 
dius the mass accretion rate M is constant (M = Mn), the ac- 
cretion timescale in the ADAF (or similar type of flow) being 
very short. Parameterizing the truncation radius as a fraction of 
the circularizat ion radius cFrank et al.-20 02) rin - fir^uc (where 
/, < 0.48, see lMenouetal.lll999h 
^ X! 0.1 one obtains 

Mn < lO-V'''^da7M0 y"', 

which because of the assumed value of q is independent of the 
primary's mass. Taking /, » 0.48 one obtains 



and assuming a mass ratio 



(7) 



Mn, <3.6xlO-^P^^;'M0y 



(8) 



which, considering the uncertainties, corresponds reasonably 
well to the observed LxiPorh) relation for quiescent SXTs 
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( iHameurv et alJ2003l) . or at least can be considered as an upper 
limit. 

Using Eq. (jS) one obtains upper limits to mass-transfer 
rates given in the last column of Table 13 They are consis- 
tent with the model B and one can see that a small amount of 
additional AML could provide an even better agreement. One 
can therefore conclude that a substantially reduced MS W AML 
would be consistent with both the stability properties and mass- 
transfer rates of SXTs. 

3.2.3. Number of transients 

Estimates of the number of SXTs in our Galaxy based on ob- 
servations range from several hundreds t o a couple tho usand 
( lvandenHeu vel"l992' IChen et al.l[l997l lRom ani"l998'). The 
number of SXTs predicted by the models is 3000 (Model A) 
to 5000 (Model B), on the upper end of the estimates based 
on observations. Assuming that their outbursts peak on av- 
erage at 0. 1 Ledd in X-rays, all of these transients should be 
seen in outburst by the RXTE ASM, regardless of their location 
in the Galaxy. The current discovery rate seems to be ^ 2-5 
per year. Matching thi s rate would impl y recurrence times of 
several hundred vears. lChen et alJ ( 119971) find that the average 
peak luminosity of transients is logL ^ 37.7 ± 0.8. The peak 
flux could be overestimated: with a peak at 10^'' erg s"' in the 
RXTE ASM bandpass, only «; half of the systems are above a 
30 mCrab threshold (appropriate for transients) after absorp- 
tion, yielding recurrence rates closer to «; 100 years. Given the 
relatively high number of transients expected, the population 
synthesis taken at face value would suggest that many tran- 
sients are sub-Eddington at maximum and/or that recurrence 
times are ^ 100 years for most sources. 

3.2.4. A population of faint LMBHB? 

A model by Menou et al. ( 1999) (MNL) may be t he answer to 
this ap parent difficulty. Following a suggestion bv lLasota et alJ 
(119961) they pointed out that if the truncation radii were sUghtly 
larger than estimated to fit observations of quiescent SXTs the 
discs would be globally stable. Indeed, stationary (M- const.) 
discs truncated at 



^in > ^c 



1.69x10' 



M 



10-iOMf 



0.375 



Ml 



0.33 



(9) 



are cold and globally stable since their accretion rates are lower 
than the critical value given by Eq. (|6} corresponding to cold 
stable equilibria. Hence, some of the systems classified by us 

as transients would instead be cold, stably accreting sources. 

1 I 1 -^ ° 

Menou e t al. ( 1999) suggested that a population of such per- 
sistent LMBHB could be present in the Galaxy. Despite their 
relatively large mass-transfer rates (M ~ 10"'" M© yr"') these 
systems would be very faint and difficult to detect so it is inter- 
esting to establish if their putative existence is compatible with 
population models. 

Figure |6l shows the distribution of stars in two population 
models with respect to stability criterion suggested by MNL. 
Here, we assume the truncation occurs at /?;„ = 0.48/?circ, 
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Figure 6. Mass transfer rates vs. orbital period for the model 
population of semidetached low-mass binaries with assumed 
black hole accretors. Upper panel - "standard" evolution model 
with disrupted magnetic braking braking (Model A). Pluses 
mark stable according to Menou et al. stability criterion for 
cold discs systems, while dots - unstable ones. Middle panel 
shows position of systems with discs stable acording to Menou 
et al. criterion in the model in which MSW does not oper- 
ate in semidetached systems with black hole accretors (model 
B). Lower panel - systems with unstable discs in model B. 
Diamonds are mass-transfer rates estimated on the base of as- 
sumed recurrence periods, like in Fig.|5l 



the radius at which overflowing mater ial from the donor star 
merges with the accretion disc (.Menou et alJ [1999 ). MNL- 
stable systems have /?,„ > Rem- The break-down over MNL- 
stable and unstable (transient) systems is given in Table |2 In 
the "standai-d" case with MSW AML (Model A) very few discs 
that were unstable vs. DlM-criterion, become stable vs. MNL- 
criterion. On the other hand, about a half of DIM unstable sys- 
tems become stable vs. the MNL criterion in the model with- 
out MSW after RLOF (Model B). This is because mass trans- 
fer rates in the GWR-only scenario are much lower than in the 
MSW-i-GWR scenario, thereby enabling the criterion for a cold 
disc to be met more easily. The number of transients for both 
models stays at ^ 2000-3000. On the other hand, truncating the 
disc at Rciic would make all transients "MNL-stable" for both 
models. It is therefore a distinct possibility that there are very 
few "classical" SXTs in the DIM sense (i.e. regularly outburst- 
ing) and that the LMBHB population divides up into hot/stable 
and cold/stable systems, the latter erupting from time to time 
due to weather variability in truncation physics or mass trans- 
fer rate. Even in this extreme case, the predicted population size 
is too small to have cold/stable systems contribute significantly 
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to the large number of faint (Lx ~ 10^' erg s ') hard X-ray 
sources seen in deep Chandra exposures toward th e Galactic 
Center JMuno et alJ2003HBelczvnski & Taaml2004 . 

3.3. Systems with very low mass ratios 

Mass transfer in systems with very low mass ratios is uncertain 
and they have been left out of the discussion until now. For val- 
ues of q ^ 0.02 the circularization radius becomes greater than 
the estimates of the outer rad ius given by Paczvhski ( 1977) and 
IPapaloizou & Pringld (119771) . Matter flowing in from the com- 
panion circularises onto unstable orbits. At q x 0.02, matter 
is added at /?circ onto orbits that can become eccentric due to 
the 3:1 resonance. At q ^ 0.005 the circularization radius ap- 
proaches the 2:1 Lindblad resonance. This can efficiently pre- 
vent mass being transferred onto the compact object (by anal- 
ogy with the gaps created by planets in protoplanetary discs 
or the inner edges of circumbinary discs) unless the timescale 
of mass transfer is such that the gap can be filled before it is 
cleared (unlikely, as the condition for gap formation is roughly 
MiIMq <: fdyn/fvis, whcrc fdyn and fvis are respectively the dy- 
namical and viscous timescales). Subsequent evolution is prob- 
ably set by the balance between GW radiation (dominating at 
low q) which drives the system inward and the torque exerted 
by the disc which drives the system outward. Regardless of this 
exact balance, since the outer radius of the accretion disc is un- 
likely to be greater than Rare, most of these systems (if accret- 
ing) are likely to be MNL stable according to the discussion in 
the preceding section. 

4. Discussion and Conclusions 

We have modeled populations of short-period semi-detached 
close binary systems containing black holes and low-mass stars 
under various assumptions about the AML mechanism. In all 
cases the obtained population can be divided into two subsets. 
One, contains short-period (< 0.1 day), very low mass-ratio 
(q < 0.02) systems; the other contains systems whose orbital 
parameters correspond to those determined for black hole bi- 
nary transients. Mass transfer in short-period systems may be 
interrupted by resonances within the primary's Roche lobe and 
their subsequent evolution and the observational consequences 
are uncertain. For the systems with mass ratios q > 0.02, the 
disk instability model provides a framework to investigate the 
observational consequences in terms of bright, persistent sys- 
tems and transients. 

Among systems with the range of periods encountered 
in SXTs (let us note that the shortest observed period of a 
LMBHB is 4. 1 hours but nothing in our model prevents such 
systems to have periods ~ 2 hours) the number of transient 
systems depends on the assumed AML mechanism. The stan- 
dard Verbunt & Zwaan (1981) formula implies the presence of 
luminous stable systems that are not observed. On the other 
hand, AML driven solely by GWR yields only transient sys- 
tems. This suggests that in short-period LMBHB the strength 
of the AML MSW mechanism is significantly reduced. Such a 
reduction is also required by the values of mass-transfer rates 
deduced from observations. Our population model shows that 



the required mass-transfer rates cannot be due to the evolved 
nature of the mass-donor. 

The number of transients in the models is in the few thou- 
sand range, at the upper end of estimates based on observa- 
tions. Better agreement would imply transients having lower 
peak fluxes and/or long recurrence times. Another possibil- 
ity is that the number of transients is reduced compared to 
these estimates. A significant fraction of transients could be 
cold and stable if the thin accretion disc is truncated in quies- 
cence. The number of such systems is highest with GWR-only 
angular momentum losses (lower mass transfer rates): with 
^trunc = OAERciic roughly half of transient SXTs then become 
cold and stable. With MSW-hGWR, the number of cold, stable 
systems is very small. This fraction can be increased to include 
all systems if the disc is truncated up to Rdw If this were the 
case, transients would be flukes (not triggered by the DIM but 
by external factors) and the regularly outbursting SXTs envi- 
sioned by the DIM would be an exception rather than the rule 
as suggested by Lasota et al. ( 199q).^ 

Re cently, IPodsiadlowski et alJ ( l2003l) and Ijustham et alJ 
( 120061) suggested that low-mass companion stars are not likely 
to provide sufficient gravitational potential to unbind during 
common-envelope phase the envelope of black-hole's massive 
progenitor. In our formulation, this is equivalent to a low value 
of the product of common-envelope efficiency and structure pa- 
rameter a A. A trial computation with aA - 0.1 resulted in a 
population with only ^ 3400 short-period binaries containing 
black holes with low-mass companions. Only the most mas- 
sive binaries avoided merger due to the fact that they lose more 
mass on the main sequence and thus have rather low-mass en- 
velopes and survive the common-envelope phase. However, 
all black holes produced in this case have mass exceeding 
^ 14 M0, substantially higher than what is observed in short- 
period LMBHB. Post-circularization period - secondary mass 
relations for these binaries are similar to these shown in Fig.|3l 
Thus, if such systems were formed, we expect evolution simi- 
lar to described in this paper, but without formation of LMBHB 
resembling observed ones. However, it is difficult to reach firm 
conclusions on this issue as both the details of the common- 
envelope phase as well as the formation of black holes are are 
very uncertain. So low-mass stars may still be able to remove 
the envelope or if indeed only the most massive binaries sur- 
vive, their black hole masses may be lower and be consistent 
with the observed masses 

We conclude that modeling the population of LMBHB, us- 
ing standard models of stellar evolution and accretion physics, 
leads to prediction of (1) a significant population of systems 
with q < 0.02 and Porb ^ 2 hours; (2) a significant number of 
persistent LMBHB unless MSW angular momentum losses are 
reduced; (3) a high number of transients unless the peak X-ray 
fluxes are low and/or the recurrence time high; (4) a possibility 
of having a large fraction, if not all LMBHB systems, actually 
being (secularly) cold and stable. 
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Appendix A: Model for the star formation history 
and stellar density distribution 

The probability of a star being born at a radius R or smaller 
depends on the integrated SFR, i.e.. 



PiRj)^ 



J^SFRiR',t)27TR'dR' 
[^"'"SFR(R',t)2nR'dR'' 



(A.l) 



where Rm^x is the maximum extent of the Galactic disc (19 
kpc). We assume that the stars do not migrate radially in 
time. The age of the Galaxy in this model is 13.5 Gyr Since 
the Boissier & Prantzos ( 1999 ) model gives the SFR projected 
onto the plane of the Galaxy, we assume a z-distribution in the 
Galactic disc 



P(z) oc sech(z/zh) 



(A.2) 



where Zh - 200 pc, neglecting dependence on age and mass. 
We have added a bulge to the model of Galactic distribution 
of stars, by d oubhng the SFR in the inn er 3 kpc of the Galaxy 
compared to Boissier & Prantzosl ( ll999l) . The total mass in this 
region at the age of 13.5 Gyr is then 2.6 10'° Mp, consistent 
with kinematic and micro-lensing results (e. g., iKlvpin et alJ 
l2002h . We distribute the stars in the bulge according to 



Pbuige = exp(-(r/0.5 kpc) ), 



(A.3) 



where r - -\Jx^ +y^ + z^. The time-dependence of SFR under 
given assumptions and compari son to often us ed exponential 
SFR may be found in Fig. 2 of JSlelemans et al.l 12004.) . 
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